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1 Introduction 

Discrete nonlinear Schrodinger (DNLS) equations can be used to model nu- 
merous phenomena in atomic, molecular, and optical physics. The general 
feature of these various settings that leads to the relevance of DNLS mod- 
els is a competition between nonlinearity, dispersion, and spatial discreteness 
(which can be periodic, quasiperiodic, or random). In three-dimensions (3D), 
the DNLS with cubic nonlinearity is written in normalized form as 

2^/,m,n — €Aui.rn,n i |^/,m,n| ^l,m.n ~\~ Vl,m,n(t*)'U>l,m,n ; (1) 

where A is the discrete Laplacian; e is a coupling constant; ui^ m ^ n is the value 
of the field at site (I, m, n); and Vi im ,n{t) is the value of the external potential 
at that site. In Eq. ([TJ) , a + sign represents the defocusing case and a — 
sign represents the focusing one. Both of these situations have been discussed 
extensively throughout this book. 

The study of DNLS equations dates back to theoretical work on biophysics 
in the early 1970s [25]. In the late 1980s, this early research motivated exten- 
sive analysis of such equations for the purpose of modeling the dynamics of 
pulses in optical waveguide arrays [17]. One decade later, experiments using 
fabricated Aluminum gallium arsenide (AlGaAs) waveguide arrays [29] stim- 
ulated a huge amount of subsequent research, including experimental investi- 
gations of phenomena such as discrete diffraction, Peierls barriers, diffraction 
management [28,60], gap solitons [54], and more [59]. As was first suggested 
theoretically in [26] and realized experimentally in [33,35,55,65], DNLS equa- 
tions also accurately predict the existence and stability properties of nonlinear 
localized waves in optically-induced lattices in photorefractive media such as 
Strontium Barium Niobate (SBN). Because of this success, research in this 
arena has exploded; structures such as dipoles [93], quadrupoles [92], multi- 
phase patterns (including soliton necklaces and stripes) [64,94], discrete vor- 
tices [34,66], and rotary solitons [89] have now been theoretically predicted 
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and experimentally obtained in lattices induced with a self-focusing nonlin- 
earity. As discussed in Ref. [84] (and references therein), self-defocusing real- 
izations have also been obtained. These allow the construction of dipole-like 
gap solitons, etc. 

DNLS equations have also been prominent in investigations of Bose- 
Einstein condensates (BECs) in optical lattice (OL) potentials, which can 
be produced by counter-propagating laser beams along one, two, or three 
directions [63]. This field has also experienced enormous growth in the last 
ten years; major experimental results that have been studied using DNLS 
equations include modulational ("dynamical") instabilities [14,81], gap soli- 
ton dynamics [27], Bloch oscillations and Landau-Zener tunnelling [4], and 
the production of period-doubled solutions [38] . 



2 Optics 

For many decades, optics has provided one of the traditional testbeds for 
investigations of nonlinear wave propagation [45]. For example, the (contin- 
uous) nonlinear Schrodinger (NLS) equation provides a dispersive envelope 
wave model for describing the electric field in optical fibers. In the presence of 
a spatially-discrete external potential (such as a periodic potential), one can 
often reduce the continuous NLS to the discrete NLS. In this section, we will 
consider some appropriate situations that arise in optical waveguide arrays 
and photorefractive crystals. Many of these optical phenomena have almost 
exact analogs in both solid state and atomic physics [79]. 



2.1 Optical Waveguide Arrays 
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Fig. 1. (Left) Experimental setup for the waveguide array experiments reported in 
Ref. [29]. (Right) Low power (diffraction) versus high power experiments. The latter 
result in discrete spatial solitons. (Images from Ref. [29].) 



An optical waveguide is a physical structure that guides electromagnetic 
waves in the optical spectrum. Early proposals in nonlinear optics suggested 
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Fig. 2. Experimental observation of a soliton (bottom panel) at the output facet 
of a waveguide array. The peak powers are 70 W (top), 320 W (center), and 500 W 
(bottom). (Image from Ref. [29].) 



that light beams can trap themselves by creating their own waveguide through 
the nonlinear Kerr effect [1] . Waveguides confine the diffraction, allowing spa- 
tial solitons to exist. In the late 1990s, Eisenberg et al. showed experimentally 
that a similar phenomenon (namely, discrete spatial solitons) can occur in a 
coupled array of identical waveguides [29]. One injects low- intensity light into 
one waveguide (or a small number of neighboring ones); this causes an ever 
increasing number of waveguides to couple as it propagates, analogously to 
what occurs in continuous media. If the light has high intensity, the Kerr effect 
changes the refractive index of the input waveguides, effectively decoupling 
them from the rest of the array. That is, certain light distributions propagate 
with a fixed spatial profile in a limited number of waveguides; see Figs. [T]and 

EI 

The standard theoretical approach used to derive a DNLS equation in 
the context of one-dimensional (ID) waveguide arrays is to decompose the 
total field (describing the envelope amplitude) into a sum of weakly coupled 
fundamental modes that are excited in each individual waveguide [17,83]. If 
one supposes that each waveguide is only coupled to its nearest neighbors, this 
approach amounts to the tight-binding approximation of solid state physics. 
For lossless waveguides with a Kerr (cubic) nonlinearity, one obtains the DNLS 

dE 

i-^+pE n + c {E n+1 +E n ^)+\\E n \ 2 E n + l i{\E n+l \ 2 +\E n _ 1 \ 2 )E n = Q, (2) 

where E n is the mode amplitude of the nth waveguide, z is the propagation 
direction, (3 is the field propagation constant of each waveguide, c is a coupling 
coefficient, and A and fi are positive constants that respectively determine the 
strengths of the self-phase and cross-phase modulations experienced by each 
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waveguide. The quantity A is proportional to the optical angular frequency 
and Kerr nonlinearity coefficient, and is inversely proportional to the effective 
area of the waveguide modes [83]. In most situations, the self-phase modula- 
tion dominates the cross-phase modulation (which arises from the nonlinear 
overlap of adjacent modes), so that /i <C A. This allows one to set ^ = and 
yields the model that Davydov employed for a-spiral protein molecules [25]. 
The transformation E n (z) = <!> n (z) exp[z(2c + j3)z] and a rescaling then gives 
a ID version of Eq. (!]). 

The impact of Ref. [29] was immediate and powerful, as numerous sub- 
sequent experiments reported very interesting phenomena. For example, this 
setting provided the first experimental demonstration of the Peierls-Nabarro 
(PN) potential in a macroscopic system [60], thereby explaining the strong 
localization observed for high-intensity light in the original experiments [29] . 
That is, the PN potential describes the energy barrier between the (stable) 
solitons that are centered on a waveguide and propagate along the waveguide 
direction and the (unstable) ones that are centered symmetrically between 
waveguides and tend to shift away from the waveguide direction. Eisenberg 
et al. [28] have also exploited diffraction management (which is analogous 
to the dispersion management ubiquitously employed in the study of tempo- 
ral solitons [53]) to produce structures with designed (reduced, canceled, or 
reversed) diffraction properties. The ability to engineer the diffraction prop- 
erties has paved the way for new possibilities (not accessible in bulk media) 
for controlling light flow. For example, using two-dimensional (2D) waveguide 
networks, discrete solitons can travel along essentially arbitrarily curves and 
be routed to any destination [18]. This may prove extremely helpful in the 
construction of photonic switching architectures. 

A waveguide array with linearly increasing effective refractive index, which 
can be induced using electro- or thermo-optical effects, has also been used 
to demonstrate Bloch oscillations (periodic recurrences) in which the initial 
distribution is recovered after one oscillation period [61]. A single- waveguide 
excitation spreads over the entire array before refocusing into the initial guide. 
More recently, Morandotti et al. investigated the interactions of discrete soli- 
tons with structural defects produced by modifying the spacing of one pair 
of waveguides in an otherwise uniform array [58]. This can be used to adjust 
the PN potential. It has also been demonstrated experimentally that even a 
binary array is sufficient to generate discrete gap solitons, which can then be 
steered via inter-band momentum exchange [62]. From a nonlinear dynam- 
ics perspective, an especially exciting result is the experimental observation 
of discrete modulational instabilities [56]. Using an AlGaAs waveguide ar- 
ray with a self-focusing Kerr nonlinearity, Meier et al. found that such an 
instability occurs when the initial spatial Bloch momentum vector is within 
the normal diffraction region of the Brillouin zone. (It is absent even at very 
high power levels in the anomolous diffraction regime.) More recent exper- 
imental observations include discrete spatial gap [15] and dark [82] solitons 
in photovoltaic lithium niobatc (LiNbOa) waveguide arrays, evidence for the 
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spontaneous formation of discrete X waves in AlGaAs waveguide arrays [50] , 
and an analog of Anderson localization (which occurs in solid state physics 
when an electron in a crystal becomes immobile in a disordered lattice) [49] . 

2.2 Photorefractive Crystals 




Fig. 3. (Left) Diagram of experimental setup for the creation of a photorefractive 
crystal lattice with electro-optic anisotropy. (Right) Typical observation of the lat- 
tice at the terminal face of the crystal. Each waveguide has a diameter of about 
7(j,m and is about 11 \im away from its nearest neighbors. (Images from Ref. [33].) 



Photorefractive crystals can be used to construct 2D periodic lattices via 
plane wave interference by employing a technique known as optical induction. 
This method, which was developed theoretically in Ref. [26] and subsequently 
demonstrated experimentally for ID discrete solitons in [35] and bright 2D 
solitons in [33], has become a very important playground for investigations 
of nonlinear waves in optics [32]. One obtains a periodic lattice in real time 
through the interference of two or more plane waves in a photosensitive mate- 
rial (see Fig. [3]). One then launches a probe beam, which experiences discrete 
diffraction (the optical equivalent of quantum tunnelling in a periodic poten- 
tial) and can form a discrete soliton provided the nonlinearity is sufficiently 
large. The model for photorefractive crystals is a continuous NLS equation 
with saturable nonlinearity [33], 

iU g + U xx + U m - 1 + ^° lu]2 U = 0, (3) 

where z is the propagation distance, (x, y) are transverse coordinates, U is 
the slowly-varying amplitude of the probe beam (normalized by the dark 
irradiance of the crystal), Eq is the applied dc field, and /; is a lattice inten- 
sity function. For a square lattice, /; = Iq sin 2 {(a; + y)/\/2} sin 2 {(a; — y)/^/2}, 
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where Iq is the lattice's peak intensity. DNLS equations have been enormously 
insightful in providing corroborations between theoretical predictions and ex- 
perimental observations (see, in particular, the investigations of discrete vor- 
tices in Refs. [16,34,66]), although they do not provide a prototypical model 
in this setting the way they do with waveguide arrays. 

For optical induction to work, it is essential that the interfering waves 
are unaffected by the nonlinearity (to ensure that the "waveguides" are as 
uniform as possible) but that the probe (soliton-forming) beam experiences a 
significant nonlinearity This can be achieved by using a photorefractive ma- 
terial with a strong electro-optic anisotropy. In such materials, coherent rays 
interfere with each other and form a spatially-varying pattern of illumination 
(because the local index of refraction is modified, via the electro-optic effect, 
by spatial variations of the light intensity). This causes ordinary polarized 
plane waves to propagate almost linearly (i.e., with practically no diffraction) 
and extraordinary polarized waves to propagate in a highly nonlinear fash- 
ion. The material of choice in the initial experiments of Ref. [33] was the 
(extremely anisotropic) SBN:75 crystal. 

The theoretical prediction and subsequent experimental demonstration of 
2D discrete optical solitons has led to the construction and analysis of entirely 
new families of discrete solitons [26,33]. As has been discussed throughout this 
book, the extra dimension allows much more intricate nonlinear dynamics to 
occur than is possible in the ID waveguides discussed above. Early experi- 
ments demonstrated novel self-trapping effects such as the excitation of odd 
and even nonlinear localized states [65]. They also showed that photorefrac- 
tive crystals can be used to produce index gratings that are more controllable 
than those in fabricated waveguide arrays. 

Various researchers have since exploited the flexibility of photorefractive 
crystals to create interesting, robust 2D structures that have the potential to 
be used as carriers and/or conduits for data transmission and processing in 
the setting of all-optical communication schemes (see Refs. [16, 32] and ref- 
erences therein). In the future, ID arrays in 2D environments might be used 
for multidimensional waveguide junctions, which has the potential to yield 
discrete soliton routing and network applications. For example, Martin et al. 
observed soliton-induced dislocations and deformations in photonic lattices 
created by partially incoherent light [55]. By exploiting the photorefractive 
nonlinearity's anisotropy, they were able to create optical structures analo- 
gous to polaron^j] from solid state physics [6] . An ever-larger array of struc- 
tures has been predicted and experimentally obtained in lattices induced with 
a self-focusing nonlinearity. For example, Yang et al. demonstrated discrete 
dipole (two-hump) [93] and quadrupole (four-hump) solitons [92] both experi- 
mentally and theoretically. They also showed that both dipole and quadrupole 
solitons are stable in a large region of parameter space when their humps are 



1 A polaron is a quasiparticle composed of a conducting electron and an induced 
polarization field that moves with the electron. 
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out of phase with each other. (The stable quadrupole solitons were square- 
shaped; adjacent humps had a phase difference of 7T, so that diagonal humps 
had the same phase.) Structures that have been observed in experiments in 
the self-defocusing case include dipole-like gap solitons [84] and gap-soliton 
vortices [44]. 

More complicated soliton structures have also been observed experimen- 
tally. For example, appropriately launching a high-order vortex beam (with, 
say, topological charge m = 4) into a photonic lattice can produce a sta- 
tionary necklace of solitons [94]. Stripes of bright [64] and gap [90] solitons 
have also been created [64], providing an interesting connection with several 
other pattern- forming systems [21]. Another very fruitful area has been the 
construction of both off-site and on-site discrete vortices [16,34,66]. Recent 
experimental observations in this direction have included self-trapping and 
charge-flipping of double-charged optical vortices (which lead to the forma- 
tion of rotating quasi- vortex solitons) [8]. Discrete rotary solitons [89] and 
discrete random-phase solitons [20] have also been observed. In fact, the re- 
sults of Ref. [20] are reminiscent of the Fermi-Pasta-Ulam (FPU) numerical 
experiments [13], as an initially homogeneous distribution in momentum space 
evolved into a steady-state multi-humped soliton power spectrum. 




Fig. 4. Experimental image of the decagonal field-intensity pattern in an optically- 
induced nonlinear photonic quasicrystal. (Image from Ref. [36].) 
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The investigation of photorefractive crystals continues to produce exper- 
imental breakthroughs, offering ever more connections to solid state physics. 
One particularly exciting experiment was the observation of dispersive shock 
waves [88] . Another fascinating result was the observation of an analog of An- 
derson localization in disordered 2D photonic lattices [78]. In this context, the 
transverse localization of light is caused by random fluctuations. Wave, de- 
fect, and phason dynamics (including discrete diffraction and discrete solitons) 
have recently been investigated experimentally in optically-induced nonlinear 
photonic quasicrystals [36,37] (see Fig. rjj, whose theoretical investigation pro- 
vides one of the outstanding challenges for DNLS models. 

3 Bose-Einstein Condensation 

At sufficiently low temperature, bosonic particles in a dilute 3D gas occupy 
the same quantum (ground) state, forming a Bose-Einstein condensate (BEC) 
[12,22,43,70]. Seventy years after they were first predicted theoretically, dilute 
(i.e., weakly interacting) BECs were finally observed experimentally in 1995 in 
vapors of rubidium and sodium [5,24]. In these experiments, atoms were loaded 
into magnetic traps and evaporatively cooled to temperatures well below a 
microkelvin. To record the properties of the BEC, the confining trap was 
then switched off, and the expanding gas was optically imaged [22]. A sharp 
peak in the velocity distribution was observed below a critical temperature 
T Cl indicating that condensation had occurred. 

If the temperature is well below T c , then considering only two-body, mean- 
field interactions, the BEC dynamics is modeled using the 3D Gross-Pitaevskii 
(GP) equation (i.e., the continuous cubic NLS equation), 



where = ^(r, t) is the condensate wave function (order parameter) nor- 
malized to the number of atoms, V(r) is the external potential, and the ef- 
fective self-interaction parameter is g — [4:Trh 2 a/m][l + 0(£ 2 )], where a is 
the two-body scattering length and ( = ^/|tf'| 2 |a| 3 is the dilute-gas param- 
eter [7,22,47]. The cubic nonlinearity arises from the nearly perfect contact 
(delta-function) interaction between particles. 

In a quasi- ID ("cigar-shaped") BEC, the transverse dimensions are about 
equal to the healing length, and the longitudinal dimension is much larger 
than the transverse ones. One can then average (H]) in the transverse plane to 
obtain the ID GP equation [22,74], 




ihu t = — [h 2 / (2m)]u xx + g\u\ 2 u + V(x)u , 



(5) 



where u, g, and V are, respectively, the rescaled ID wave function, interac- 
tion parameter, and external trapping potential. The interatomic interactions 
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Fig. 5. Diagram of the creation of a ID optical lattice potential using (a) counter- 
propagating laser beams and (b) beams intersecting at angle 6. The quantities 
and k' L denote the wave vectors of the beams. The lattice period is given by the 
distance d between consecutive maxima of light intensity in the interference pattern. 
(Images from Ref. [63].) 



in BECs are determined by the sign of g: they are repulsive (producing a 
defocusing nonlinearity) when g > and attractive (producing a focusing 
nonlinearity) when g < 0. 

BECs can be loaded into optical lattice (OL) potentials (or superlattices, 
which are small-scale lattices subjected to a large-scale modulation), which 
are created experimentally as interference patterns of laser beams. Consider 
two identical laser beams with parallel polarization and equal peak intensi- 
ties, and counterpropagate them as in Fig. [5^i so that their cross sections 
overlap completely. The two beams create an interference pattern with period 
d = Al/2 (half of the optical wavelength) equal to the distance between con- 
secutive maxima of the resulting light intensity. The potential experienced by 
atoms in the BEC is then [63] 
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Fig. 6. Experimental demonstration of a bright gap soliton (the shading shows 
the atomic density). A small, stable peak forms after about 25 ms. (Image from 
Ref. [27].) 

V{x) = V cos 2 (7rx/d), (6) 

where Vo is the lattice depth. See Ref. [63] for numerous additional details. 

BECs were first successfully placed in OLs in 1998 [4], and numerous labs 
worldwide now have the capability to do so. Experimental and theoretical in- 
vestigations of BECs in OLs (and related potentials) have developed into one 
of the most important subdisciplines of BEC investigations [63,70,71]. We 
focus here on results that can be modeled using a DNLS framework; see the 
reviews [31,63] for discussions of and references to myriad other outstanding 
experiments. The first big experimental result was the observation of Bloch 
oscillations in a repulsive BEC by Anderson and Kasevich [4]. They used a 
trapping potential with both an OL and a linear (gravitational) component 
to create a sloping periodic ("washboard") potential. When the slope was 
small, wave packets remained confined in a single band and oscillated coher- 
ently at the Bloch frequency. This effect is closely related to the ac Josephson 
effect in superconducting electronic systems. For larger slopes, wave packets 
were able to escape their original band and transition to higher states. Trom- 
bettoni and Smerzi analyzed the results of these experiments using a DNLS 
equation that they derived from the CP equation with the appropriate (OL 
plus gravitational) potential using the tight-binding approximation valid for 
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moderate-amplitude potentials in which Bloch waves are strongly localized in 
potential wells [85]. In this paper, they also showed that discrete breathers 
(specifically, bright gap solitons) can exist in BECs with repulsive potentials 
(i.e., defocusing nonlinearities) . 

Other fundamental experimental work on BECs in OLs has concerned 
superfluid properties. In 2001, Burger et al. treated this setting as a homo- 
geneous superfluid with density-dependent critical velocity [11]. Cataloiotti 
et al. then built on this research to examine a classical transition between 
superfluid and Mott insulator behavor in BECs loaded into an OL superim- 
posed on a harmonic potential [14]. (The better-known quantum transition 
was first shown in Ref. [41].) The BEC exhibits coherent oscillations in the 
"superfluid" regime and localization in the harmonic trap in the "insulator" 
regime, in which each site has many atoms of its own and is effectively its 
own BEC. The transition from superfluidity to Mott insulation occurs when 
the condensate wave packet's initial displacement is larger than some critical 
value or, equivalently, when the velocity of its center of mass is larger than 
a critical velocity that depends on the tunnelling rate between adjacent OL 
sites. These experiments confirmed the predictions of Ref. [81], which used a 
DNLS approach to predict the onset of this superfluid-insulator transition via 
a discrete modulational ( "dynamical" ) instability and to derive an analytical 
expression for the critical velocity at which it occurs. 

Subsequent theoretical work with DNLS equations predicted that modu- 
lational instabilities could lead to "period-doubled" solutions in which the 
BEC wavcfunction's periodicity is twice that of the underlying OL [52]. 
(Period-doubled wavefunctions were simultaneously constructed using a CP 
approach [72].) The modulational instability mechanism was exploited experi- 
mentally the next year to construct these solutions by parametrically exciting 
a BEC via periodic translations (shaking) of the OL potential [38] . Parametric 
excitation of BECs promises to lead to many more interesting insights in the 
future. 

By balancing the spatial periodicity of the OL with the nonlinearity in the 
DNLS, one can also construct intrinsic localized modes (discrete breathers) 
known as bright gap solitons, which resemble those supported by Bragg grat- 
ings in nonlinear optical systems. In BECs, such breathers have been predicted 
in two situations: 

1. The small-amplitude limit in which the value of chemical potential is close 
to forbidden zones ( "gaps" ) of the underlying linear Schrodingcr equation 
with a periodic potential [48]. 

2. In the tight-binding approximation, for which the continuous NLS equa- 
tion with a periodic potential can be reduced to the DNLS equation [81]. 
(As mentioned earlier, this corresponds to the standard manifestation of 
the DNLS in waveguide arrays.) 

Recent experiments [27] have confirmed the first prediction (see Fig. [6]). 
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Another important development was the experimental construction of 2D 
and 3D OLs [63], which as in optical systems leads to much more intricate 
nonlinear localized structures. One can obtain higher-dimensional OL poten- 
tials by using additional pairs of laser beams. The simplest way to do this is 
to have pairs of countcrpropagating laser beams along each of two or three 
mutually orthogonal axes. The interference pattern obtained with this many 
laser beams depends sensitively on their polarizations, relative phases, and 
orientations. This allows experimentalists to construct a large variety of OL 
geometries in 2D and 3D. In 2001, Greiner et al. showed that BECs can be 
efficiently transferred into 2D lattice potentials by adiabatically increasing the 
depth of the lattice [39,40]. They confined atoms to an array of narrow poten- 
tial tubes, each of which was filled with a ID quantum gas. Around the same 
time, Burger et al. confined quasi-2D BECs into the lattice sites of a ID OL 
potential [10]. By adding more laser beams and/or controlling their polariza- 
tions and relative phases, experimentalists can in principle create even more 
complicated potentials (such as quasiperiodic or Kagome lattices) [75]. Very 
recently, there has been also been a great deal of theoretical and experimental 
interest in rotating optical lattices [86], which can likely be modeled using 
an appropriate DNLS framework to study interesting vortex dynamics. Such 
experiments and (more generally) investigations of BECs with effective fields 
obtained by this and other [42] means promise to yield considerable insights 
into quantum hall physics. 

In addition to the fascinating insights into basic physics discussed above, 
two other major consequences of investigations of BECs in OLs and related 
potentials have been to bring quantum computation one (small) step closer 
to reality and to help bridge the gap between condensed matter physics and 
atomic/molecular physics. One of the key proposed systems for constructing 
a quantum computer is a BEC in optical lattice and related potentials [19,87]. 
This was the motivation for the experimental implementation of optical su- 
perlattice potentials [69] and its 2D egg-carton descendent, which consists of 
an optical lattice potential in one cardinal direction and a double-well poten- 
tial in another [3]. This has led very recently to the experimental realization 
of a two-qubit quantum gate [2] . The second front has been advanced exper- 
imentally in the OL context by investigations of Fermi condensates in OL 
potentials [46,57], Bose- Fermi mixtures in OL potentials [68], and more. Fi- 
nally, in parallel with the recent insights in optics discussed above, Anderson 
localization has been observed recently both for a BEC placed in a ID waveg- 
uide with controlled disorder [9] and for a BEC in a ID quasiperiodic OL [73]. 

4 Summary and Outlook 

In this review, I have discussed applications in nonlinear optics and Bose- 
Einstcin condensation in which DNLS equations have been used to explain 
fundamental and striking experimental results. In optics, DNLS equations pro- 
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vide a prototypical model for the dynamics of discrete solitons in waveguide 
arrays. The same is true for BECs in optical lattice (and related) potentials. 
DNLS equations have also been used successfully to predict robust experi- 
mental features in photorefractive crystals, although they do not provide a 
prototypical model in this setting. 

DNLS equations arise in a number of other contexts as envelope models 
for several types of nonlinear lattice equations (such as ones of Klein-Gordon 
type). Related experiments have revealed the existence of intrinsic localized 
modes in these systems [76]. Relevant settings include quasi-lD antiferro- 
magnets [91], micromechanical oscillator arrays [76], and electric transmission 
lines [77]. 

DNLS equations have also been used in a variety of other settings to 
make interesting predictions that have not yet been verified experimentally. 
For example, in composite metamaterials, Shadrivov et al. have analyzed the 
modulational instability of different nonlinear states and demonstrated that 
nonlinear metamaterials support the propagation of domain walls (kinks) that 
connect regions of positive and negative magnetization [80] . Very recently, dis- 
sipative discrete breathers were constructed in a model of rf superconducting 
quantum interference device (SQUID) arrays [51]. (Similar discrete breathers 
have also recently been studied theoretically in both ID and 2D in the setting 
of metamaterials [30].) This model is reminiscent of a DNLS with dissipation, 
except that the nonlinearity was sinusoidal rather than cubic. A bit farther 
afield, ID chains of granular materials (sometimes called "phononic crystals" ) 
have been given increasing attention from both experimentalists and theorists 
in recent years [23,67]. When given an initial precompression, they can ex- 
hibit optical modes that are expected to be describable as gap solitons in a 
nonlinear lattice model reminiscent of FPU chains. 

In conclusion, DNLS equations and related models have been incredibly 
successful in the description of numerous experiments in nonlinear optics and 
Bose-Einstein condensation. They also show considerable promise in a number 
of other settings, and related nonlinear lattice models are also pervasive in a 
huge number of applications. To borrow a phrase from the defunct rock band 
Timbuk3, the future's so bright that we've got to wear shades. 
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